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SUMMARY 


PROBLEM 

'  ^The  objective  of  this  work  was  to  develop  a  possible  way  of  source  ranging  by  a 
single  array  in  a  multipath  ocean  and  to  indicate  its  limitations. 

RESULTS 

It  was  shown  that  mean-square  sound-level  fluctuations  for  a  single  sensor  in  a 
random  multipath  environment  are  directly  dependent  on  the  source  range.  When  this  was 
applied  to  an  array  of  sensors,  the  result  showed  that  a  single  horizontal  line  array  might  track 
a  moving  source  in  a  multipath  ocean  environment.  This  result  would  facilitate  the  inter¬ 
array  communication  problem  for  coherent  multi-array  processing  (CMAP)  tracking  of  sub- 
jiiarines.^ln  Section  IV  of  Ref  1  this  single-array  source  tracking  method  was  tested  with 
promising  results  for  the  Bearing  Stake  data. 

RECOMMENDATIONS 

It  is  recommended  that  this  tracking  method  be  further  tested  with  even  more 
suitable  data. 


1 .  Reference  available  to  qualified  requestors. 
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NOTATION 


amplitude  of  acoustic  signal  F 
phase  of  acoustic  signal  F 
sound-level  of  Eq  ( 1 5) 
acoustic  signal  on  path  m  at  sensor  j 
source  bearing 
wave  number 

the  mean-square  fluctuation  of  the  refractive  index 

the  integral  scale  for  the  correlation  function  of  the  refractive  index 

amplitude  for  a  deterministic  medium 

phase  for  a  deterministic  medium 

mean-square  sound-level  fluctuation  for  sensor  j 

mean-square  sound-level  fluctuation  for  the  array 
statistically  independent  noise  power 


I.  INTRODUCTION 


References  2,  3,  and  4  derived  and  verified  useful  relations  for  the  mean-square 
sound-level  fluctuation  and  the  mean-square  phase  fluctuation  for  single  paths  to  a  single 
sensor  in  the  ocean.  In  this  report,  these  relations  are  extended  to  give  the  mean-square 
sound-level  fluctuation  for  a  single  sensor  in  a  random  multipath  environment.  This  rela¬ 
tion  could  permit  source  ranging  by  a  single  sensor,  and  when  applied  to  an  array  of 
sensors,  the  result  shows  that  a  single  horizontal  line  array  may  track  a  moving  source  in 
a  multipath  ocean  environment. 


2.  V.  I.  Tatarski,  Wave  Propagation  in  a  Turbulent  Medium,  McGraw-Hill,  195 1 . 

3.  R.  F.  Shvaehko,  "Sound  Fluctuations  in  the  Upper  Layer  of  the  Ocean  and  Their  Relation  to  the  Random 
Inhomogeneities  of  the  Medium,"  .Vor.  Plivs.-Acoiist.  _9,  280-282  (1964). 

4.  R.  F.  Shvaehko,  "Sound  Fluctuations  and  Random  Inhomogeneities  in  the  Ocean."  Sow  Phvs.-Acoust. 

_L3.  93  -97  ( 1967).  ; 
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11.  STATISTICAL  CHARACTERISTICS  OF  MULTIPATH  SIGNALS 


The  multipath  signals  received  by  a  horizontal  line  array  can  be  represented  in  a 
multipath  ocean  environment  as  follows.  First  assume  that  an  ocean  acoustic  signal  F  re¬ 
ceived  at  sensor  point  x  and  time  t  can  be  represented  by 

F  =  F(x,t)  =  Ae#  ,  (1) 

where  A  is  the  amplitude  and  <j>  is  the  phase.  Next,  assume  that  F  is  composed  of  a  linear 
superposition  of  multipath  signals  Fnr  m  =  1,  2,  3 . M  (M  is  the  total  number  of  multi¬ 

path  signals  received)  and  that  each  multipath  signal  received  can  be  represented  by 


Fm  = 


Ame 


¥ 


m 


(2) 


The  j**1  sensor  of  a  horizontal  line  array  with  sensors  at  Xj,  j  =  1 , 2, 3, .  . . ,  J  (J  is  the  total 
number  of  array  sensors)  receives  a  signal  represented  by 


id)-. 

F-  =  A  e  J 
rJ  rtJ 


M 


m=l 


M 

I 

m=l 


=  V  F  =  y  A  e^mj 

Z,  mj  Z. 


(3) 


(4) 


Note  that  Eq  (4)  represents  an  irreversible  linear  superposition  at  the  array  in  the 
sense  that  the  multipath  sig'nal  components  Fmj.  not  to  mention  their  individual  amplitudes 
Amj  and  phases  0mj,  cannot  be  resolved  by  any  processing  techniques  on  the  array  signal  Fj 
For  example,  take 

M  M 

lFj|2  =  Af  =  X  X  AmjAnjcos(,»mj- V  •  (5) 

m=l  n=l 


and  note  that  not  only  are  the  component  amplitudes  mixed  together,  but  also  that  the  com¬ 
ponent  phases  remain  in  a  manner  determined  by  the  algorithm  used  for  determining 
Likewise, 


cot  <t> j 


Re  Fj 

ImTj 


M 

X  Amj  cos  ^mj 
m=l 
M 

£  Amj  s‘n  ^mj 
m=l 


(6) 


illustrates  that  the  component  amplitudes  remain  in  0j  in  a  manner  that  is  determined  by  the 
algorithm  employed  for  determining  <t> j . 
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For  later  convenience,  define  the  source  bearing  4>0  as  the  appropriate  peak  B(4>0) 
of  the  phase-only  beam  former 

J  J 

B(<t>)  =  -j  ^  ^  <cos  -0e)  -  kQ  (dj  -dg)cos4>]  >  ,  (7) 

J"  j=l  fi=l 

where  kQ  is  the  wave  number,  dj  is  the  distance  from  the  first  sensor  to  the  sensor  (dj  =  0) 
and  the  operator  <•>  represents  a  time  average.  The  source  bearing  <h0  is  defined  as  zero  for  a 
forward  endfire  arrival. 

Now  consider  the  statistical  characteristics  of  multipath  signals  received  by  horizontal 
line  arrays  in  the  ocean.  Assume  the  length  of  the  array  dj  is  much  less  than  the  range  R  to 
the  array,  ie,  dj  <<  R.  Since  the  array  is  horizontal,  the  m1*1  paths  to  sensors  j  and  C,  j  =#  !<! , 
are  approximately  parallel  in  the  vicinity  of  the  array  (since  the  paths  are  of  the  same 
propagation  type)  so  that  Fmj  and  Fmg,  j  =£  £,  are  statistically  dependent  in  the  sense  that 


Amj  is  statistically  dependent  on  Aj^g.j^C  (8) 

j  is  statistically  dependent  on  j  ¥=  2  (9) 

Define  the  residual  phase  dm j  so  that 

*mj  H  %  +  kodjcos<lV 

Now  assume  that 

0m j  and  0nj  are  statistically  independent,  m  =A  n  (11) 

Amj  and  0ng  are  statistically  independent,  m  =£  n  (12) 

Amj  and  Anj  are  statistically  independent,  m  ¥=  n  .  (13) 


7 


III.  TATARSKI’S  RELATIONS  FOR  SOUND-LEVEL  AND 
PHASE  FLUCTUATIONS 


Modify  the  development  starting  on  page  124  of  Ref  2  as  follows.  Represent  the 
pressure  wave  p  solution  of  the  stochastic  Helmholtz  equation  as 


p  =  F  =  A  e'^j  =  eXmJ  +  ‘^mi 
p  rmj  Amj e  e 

where  the  sound-level  is  defined  as 

X  =  fcn  A 

and  use  the  accurate  and  practical  approximation  for  the  time  average  of  x, 


(14) 


(15) 


mj 


Xmj  Amj  A0  •  * 

where  A0  represents  the  amplitude  for  a  deterministic  medium  (with  the  subscripts  m  and  j 
suppressed).  Define  the  sound-level  fluctuation  by 

xmj  =  xmj  ~  xmj  “  (Amj/A-0)  • 

From  Ref  1  through  3, 

(xmj)-  ~  {®n  (Amj/'A0)1‘  =  p2  Ln  kQ2  R.  V/XR_»  L0  .  (18) 

where  p-  is  the  mean-square  fluctuation  of  the  refractive  index,  Ln  is  the  integral  scale  for 
the  correlation  function  of  the  refractive  index  fluctuation.  X  is  the  wavelength,  and  L0  is 
the  longest  fluctuation  length  scale  of  interest.  Reference  3  found 

P2  Ln  ~  4.2  •  I0~9  m  ( 14) 

in  the  Atlantic  Ocean.  An  accurate  and  practical  approximation  for  the  time  average  of 
Omj  is 


^mj  "  *0  "  ko  dj  cos  *o  •  (20) 

where  0O  represents  the  phase  tor  a  deterministic  medium  (with  the  subscripts  m  and  j  sup¬ 
pressed)  and  k0  dj  cos  <l>0  is  the  phase  delay  from  the  first  to  the  j1*1  sensor.  Equation  ( 20) 
is  an  appropriate  relation  for  array  studies.  Define  the  phase  fluctuation  as 

^  mj  =  ^mj  ~  ^mj  *  - '  * 


^mj  *co^jcos^’o  ^mj 


(22) 
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Equation  (9.35)  of  Ref  2  gives 


6 “  • 
mj 


(0mj  ^o>"  ^o-  Lo  ' 


(23) 


q 


IV.  DEVELOPMENT  OF  THE  MEAN-SQUARE  SOUND-LEVEL 
FLUCTUATION  FOR  A  MULTIPATH  MEDIUM 


The  mean-square  sound  fluctuation  relation,  ( Xj)“ -  for  a  single  array  sensor  in  a  ran¬ 
dom  multipath  medium  can  be  developed  in  a  straightforward  manner  by  using  the  expression 


=  Kn  AJ 


M  M 

=  2  Y  Kn  Amj  + 

n~  1 


M  M 


Kn  Anj  + 


^  tin  cos  <0mj-0nj>- 


m=  I 


m=l  n=l 


which  comes  from  F.q  (5)  and  (15).  Therefore, 


(Xj>“  =  (Xj  -  Xj)“ 


(Sin  Aj)“  -  ((in  Aj)- 


M  M 


V 

z- 

y  ( 

«"Amj 

m=l 

n=  1 

M 

M 

M  , 

V 

z- 

V 

2b 

m'= 

1  m  =  1 

n=  i 

M 

M 

M 

V 

V 

V 

m= 

1  11=1 

m'=  1  i 

nn.i  m 


Sin  Am'j  Sin  cos(0ny  -0nj)  -  Cn  AnVj  C.i  cos(0mj  -  <?nj) 


( Sin  lcos(0mj  —  >1  Cn  ( cos -0n'j>l 

-  tin  [cos  <0mj -0nj»l  Sin  Icos  (0m-j  -  0n-j )]  )  . 


( 25 ) 


The  last  term  in  Lq  (25  I  is  approached  as  follows.  Applying 


£n  cos  (0mj  -0nj) 


n j  u  n j ' 


{  ~  _  i(P  .  +  p  \ 

I  2  \  mj  nj  / 

I’0 


£n  [cos  (0mj-  -0nj)l  £n  [cos  (tf’m'j  ~  4>n'j)l  =  6n  l  cos  (&mj  -  0nj)l  £n  Icos  (0m'j  -  0n-j)] 

(1  ^  T*  ^  ^  t  T"  \  ,  , 

—  -r  6~  .  6~ +  d~  .  6  +  6~.  9~ +  0“.  0“.  ) .  m  =£n,  m  =£ n  , 

4  mj  mj  mj  nj  nj  mj  nj  nj/ 

=  j  (29 

1=0  ,  m  =  n  or  m'  =  n' , 

via  Eq  ( 1  1). 

Now  consider  the  last  term  in  Eq  ( 25).  as  follows: 

JV1  M  M 

(  «n  [cos  (0mj  -  0nj))  tn  Icos(0m7-0n7)l 


m=  I  n=  1  m  =  1  n : 


-  £n  [cos  (<&mj  -0nj)I  £n  [cos  (0m’,  -0n'j'l  ) 


mj  rnj 1 


Mr  1 

-  4  y  K-iw.y 

2  mj  \  mj/ 

m=l  L 


=  -VM  3(0-  )  -  <0"  ■) 

2  mj  mj 


=  5M  (0“ .)"  .  (32) 

mj 

=•  5M  (/u~  Ln  k”  R)'  .  v/XR  »  L0  .  (33) 

Equation  (30)  follows  from  Eq  (II).  (28).  and  (2l>).  The  central  limit  theorem  yields  Eq  (31). 
and  Eq  (23)  gives  Eq  (33). 

Equations  ( I  2)  and  ( 26)  give 


£n  Am-j  £n  (cos  (0mJ  -0nj)l  =  £n  Am-j  £n  [cos  (0~  -  0nj)l 


=  £n  Am-j  Sn  [cos  (0)T)j  -0n-j)[ 


so  that  the  second  term  in  Eq  (35)  goes  to  zero.  Applying  Eq  (13)  yields 


^mj  ^nj  ^mj  ^nj 

which  simplifies  the  first  term  in  Eq  (29). 

Therefore,  applying  Eq  (15),  (17),  (18),  (23),  (30),  (34),  and  (35)  to  Eq  (25)  pro¬ 
duces  the  mean-square  sound-level  fluctuation  for  a  single  array  sensor  in  a  random  multi- 
path  medium: 


M 

“  “ 

M 

-  “ 

(x  f-  =  4  ! 

n=l 

(CnAmj)2  -  (Cn  Amj)2 

+  !  I 

m=l 

^-(d2.) 2 
nj  mj 

-  I 


M 

I 

m=l 


4(x'  •)-  +  5  (O'.)' 
VA-mj  mj' 


M 


Ln  kQ-  R  (4  +  5  Ln  kQ-  R) 
—  4M  Ln  kQ“  R .  /XR~»L0, 


,  v/\R  »  L0  . 


(36) 


(37) 


(38) 


(39) 


when 

1  »  pi2  Ln  kQ~  R  .  (40) 

References  1 .  2.  and  3  give  other  relations  that  are  valid  when  the  inequality  ,/XR  »  LQ 
is  not  satisfied.  In  these  cases,  the  analysis  that  led  to  Eq  (39)  can  be  extended.  From  Eq 
( 19)  and  (38).  it  is  apparent  that  there  is  a  iarge  range  of  application  for  Eq  (39)  for  low  fre¬ 
quencies.  Shvachko  (Ref  3  and  4)  has  shown  that  stochastic  parameters  like  Ln  are  stable 
in  the  ocean  (see  also  Ref  5).  Therefore,  when  the  sound  propagation  conditions  are  well 
understood  (so  that  M  and  pi-  Ln  are  known),  relations  like  Eq  (39)  can  be  used  to  determine 
the  source  range  R  from  measurements  of  (Xj  )*- .  Thus,  passive  ranging  by  a  single  sensor  is 
possible  in  a  random  multipath  medium.  Equations (38)  and  (39)  apply  when  there  exists 
only  one  dominant  source,  since  they  depend  on  only  one  range. 


5.  J.  A.  Neubert.  "Experimental  Agreement  of  Stochastic  Ray-Theory  Relations,"  .1  Acnusi.  Snr.  Am. 
62.  326-334(1977). 


V.  CONCLUSIONS 


Equation  (39)  can  be  applied  to  a  horizontal  line  array  so  that  the  source  bearing  4>0 
is  available  as  well  as  the  range  R.  This  could  permit  source  tracking  by  a  single  horizontal 
line  array  in  a  multipath  ocean  environment.  The  total  array  gives 

_  J 

(x')“  =  ^  (Xj)2  +  Nq  (41) 

j=l 

-  4JMM2Lnk02R  +  N0,  ,/XR  »  L0  (42) 

for  F.q  (39).  In  Eq  (41)  and  (42)  NQ  has  been  included  to  show  that  statistically  independent 
noise  N0  only  results  in  a  range-independent  offset  to  the  basic  source-ranging  relation  when 
one  dominant  source  is  embedded  in  a  random  noise  field.  Equation  (42)  shows  that  the 
mean-square  sound-level  fluctuation  level  increases  by  the  number  J  of  sensors  when  their 
outputs  are  summed.  Hence  a  horizontal  line  array  may  determine  the  source  range  R  even 
when  the  fluctuation  level  of  a  single  sensor  is  insufficient. 

Although  Section  IV  of  Ref  1  showed  promising  results  for  this  array  source  tracking 
method  when  Bearing  Stake  data  were  used,  there  are  some  definite  limitations  that  must  be 
considered. 

1.  Although  the  lower  limit  on  range  should  present  no  practical  difficulties,  upper 
range  limits  such  as  Eq  (40)  will.  In  fact,  the  fluctuations  will  always  have  a  saturation 
limit  as  the  range  increases. 

2.  The  method  is  restricted  to  a  relatively  simple  multipath  environment.  Bottom 
loss  and  refractive  effects  can  invalidate  relations  like  Eq  (42)  since  these  latter  cannot  treat 
a  variable  number  M  of  arrivals.  By  a  different  method.  Ref  0  treats  these  and  related  pro¬ 
blems. 

3.  The  theory  behind  this  method  (Sections  III  and  IV)  can  only  treat  the  first  term 
in  Eq  (41 )  and  then  only  for  the  presence  of  a  single  dominant  source  (since  this  method 
cannot  beam  form).  References  I  and  6  show  that  N0  of  Eq  (41 )  plays  a  dominant  role. 
Although  N0  can  be  measured,  it  cannot  be  predicted,  so  unmeasured  changes  in  N0  militate 
against  the  use  of  this  source-ranging  method. 

4.  As  seen  in  Ref  1  and  (>.  this  source-ranging  scheme  requires  considerable  statistical 
stability  and.  therefore,  may  not  be  feasible  unless  suitable  processing  techniques  are  care¬ 
fully  applied. 


(>.  Reference  available  to  qualified  requestors. 
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